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An imaging system, an X-ray imaging apparatus, a method and a computer program arranged 
for reducing an artifact in a three-dimensional reconstructed volume 



The invention relates to an imaging system arranged to reduce an artifact in a 
three-dimensional reconstructed volume comprising a plurality of planar images. 

The invention further relates to an X-ray imaging apparatus comprising an 
imaging system arranged to reduce an artifact in a three-dimensional reconstructed volume. 
5 Hie invention still further relates to a method for reducing an artifact in a 

three-dimensional reconstructed volume comprising a plurality of planar images. 

The invention still further relates to a computer program arranged to carry-out 
the steps of the above method. 



10 

An embodiment of an imaging system arranged to reduce an artifact in an 
image is known from EP 0 562 657 A2. The known imaging system is used in the field of X- 
ray diagnostics, where a patient to be examined is positioned between an X-ray source and an 
X-ray detector arranged to substantially absorb transmitted through the patient X-rays. The 

1 5 X-ray detector comprises an image intensifier arranged to amplify a measured signal. In order 
to obtain an angulated projection, the X-ray source and the X-ray detector together with the 
image intensifier are rotated about a rotational axis. The known imaging system is arranged 
to reduce the image artifact occurring at the edges of the image, known as vignetting, said 
artifact being caused by a local variation in a gain factor for pixels in an acquisition system. 

20 In the known imaging system compensating for vignetting is achieved by multiplying pixel- 
values of the vignetted image by a gain characteristic yielding from a calibration 
measurement to obtain an unperturbed image. The gain characteristic is obtdned from an 
image of a calibration object with a known homogeneous brightness distribution. 



25 

It is a disadvantage of the known imaging system that a residual image artifact 
persists in the image corrected for vignetting. In particular, rings appear in slices 
reconstructed from a three-dimensional data, said slices being spatially oriented 
perpendicular to the rotational axis around which the image intensifier is rotated during a data 
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acquisition for purposes of the three-dimensional reconstruction. This ring-shaped artifact is 
caused by a difference in a local response of the image intensifier of the imaging system of 
the X-ray apparatus. 

5 

It is an object of the invention to provide the imaging system where the ring- 
shaped artifact is substantially reduced. 

The imaging system according to the invention comprises image artifact 
reduction means arranged to subsequently process said planar images with a &st corrective 

10 image for eliminating a first source of structured noise in said images and a second corrective 
image for eliminating a second source of structured noise in said images. 

The technical measure of the invention is based on the insight that gain 
differences for the image intensifier are mainly caused by inhomogeneities in the input and 
output screens of the image intensifier. Such response differences are further referred to as 

1 5 structured noise. The image intensifier suffers from image deformation due to pincushion 
distortion and due to external magnetic fields, this deformation being dependent on the 
orientation of the image intensifier. This deformation is referred to as warping. A warped 
image can be corrected in order to remove warping, said correction being carried out by 
means of an application of a per se known unwarping function; this being referred to as 

20 unwarping. During a three-dimensional scan, the relation between an element on a projection 
screen , for example a CCD camera of the image intensifier, and a path of an associated X- 
ray changes due to orientation-dependent image deformation. In order to substantially reduce 
the image artifact caused by local gain differences in both the input and output screens of the 
image intensifier a double-stage calibration is performed, by means of the first corrective 

25 image and the second corrective image, which independently account for structured noise of 
the input screen and the output screen of the image intensifier. 

In an embodunent of the imaging system according to the invention, the first 
corrective image comprises a first gain correction data, the second corrective image 
comprises a second gain correction data. 

30 It is advantageous to use gain correction data for constructing the respective 

corrective images as the structured noise of the input screen and the output screen of the 
image intensifier is reflected in local gain differences of respective screens. The first 
corrective image accounts for the structured noise at the output screen of the unage intensifier 
and the second corrective image accounts for the structured noise at the input screen of tiie 
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image intensifier. The sources of the structured noise together with corresponding suitable 
corrective images will be explained in detail with reference to Figures 3. 

In a further embodiment of the imaging system the first gain correction data 
comprises a result of an averaging of a plurality of raw images of a gain calibration scan; the 

5 second gain correction data comprises a result of an averaging of a plurality of raw images of 
the gain calibration scan after them being processed by means of the furst corrective image 
and an unwarping function. It is found to be advantageous to construct the first gain 
correction data from the raw images of the gain calibration scan, which is acquired prior to 
acquisition of, for example, a patient data scan. The first gain correction is applied to the raw 

10 images of the gain calibration scan yielding a set of gain-corrected images of the gain 

calibration scan, where the high jfrequency part of the structured noise of the input screen is 
removed. These images are then processed with a suitable unwarping function after which 
these unages are averaged yielding the second gain correction data corresponding to the 
structured noise of the input screen of the image intensifier. This embodiment is 

15 advantageous as it presents a simple possibility to carry-out a double-stage gain correction 
correcting for the structured noise of the mput and tiie output screens of the image intensifier 
separately. 

A method for reducing an artifact in a three-dimensional reconstructed volume 
comprising a plurality of planar images, according to the invention comprises the steps of: 
20 - processing the image using a first corrective image to eliminate a first source of structured 
noise in the image yielding a first-stage corrected image; 

- processing the first-stage corrected image with a second corrective image to eliminate a 
second source of structured noise in tiie image. 

25 

These and other aspects of the invention will be discussed with reference to 

figures. 

Figure la shows a schematic view of an X-ray imaging apparatus as known in 

the art. 

30 Figure lb shows a schematic view of an embodiment of an imaging system of 

the image intensifier as is used in the X-ray apparatus known in the art. 

Figure 2a shows an image artifact in a slice perpendicular to the rotational axis 
of the X-ray system of Figure la post first gain correction as known in the art. 
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Figure 2b shows an image artifact in a slice perpendicular to the rotational axis 
of the X-ray system of Figure la post second gain correction as known in the art. 

Figure 3 shows a schematic view of a block diagram of an embodiment of the 
imaging system according to the invention. 
5 Figure 4 shows schematically a flowchart of a computer program according to 

the invention. 

Figure 5 shows a resulting image in a slice perpendicular to the rotational axis 
of the X-ray system of Figure 3 post double-stage gain correction, according to the invention. 

10 

Figure la shows a schematic view of an embodiment of an X-ray imaging 
apparatus as is known in the art. The X-ray apparatus 1 comprises a C-arm lb on which an 
X-ray source Ic and the X-ray imaging system Id are mounted. The C-arm lb is rotatably 
arranged on a stand la. During an acquisition of a three-dimensional scan, the C-arm lb is 

15 rotated around it's rotational axis le. The X-rays emanating from the X-ray source Ic are 
intercepted by the patient (not shown), who is located in an imagmg volume between the X- 
ray source Ic and the X-ray detector Id. The X-rays If emanating from the X-ray source, 
which are subsequently transmitted through the patient are absorbed by the X-ray unaging 
system Id, which is discussed in detail with reference to Figure lb. 

20 Figure 1 b shows a schematic view of an embodiment of an imaging system 1 5 

of the image intensifier as is used in X-ray systems known in the art. For medical diagnostics 
purposes it is conventional to examine a patient in question by means of X-rays. For this 
purpose the patient 3 is located in an examination space (not shown) of the X-ray apparatus 1 
and a set of images with different projection angles is acquired for purposes of a three- 

25 dimensional reconstruction. The transmitted X-rays 2 intercept the input screen of the image 
intensifier 4, The electronic signal is amplified in the image intensifier and is made available 
to the optics 7 at the output screen 6 of the image intensifier. The light beam is converged at 
an input of a suitable projection means 8, for example a CCD camera or a CMOS-detector. 
An output signal from the CCD camera 8 is made available for fiirther processing at the 

30 output 9. There are several sources of noise in the X-ray system. First, the X-ray beam 2 

comprises an X-ray noise 2a, which is mainly caused by a quantum noise of the X-ray source 
and contribution of a scatter radiation. The image intensifier has two main sources of noise: a 
structured noise 4a of the input screen and a structured noise 6a of the output screen of the 
image intensifier. The optics 7 contributes to the overall noise by its structure noise and a 
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shot noise 7a. Additionally the CCD camera 8 has a readout noise component and an 
acquisition noise component. It is a common practice to calibrate the CCD camera for gain 
and of&et, therefore the remaining m^n noise components in the output signal at output 9 
comprise the structured noise at the input and output screens of the image intensifier. The 

5 contribution of these noise sources to the artifact in a three-dimensional reconstructed volume 
is shown in detail with reference to Figures 2a and 2b. 

Figure 2a shows an image artifact post a per se known gain correction for a 
slice II corresponding to a plane perpendicular to the rotational axis le of the image 
intensifier as is set forth with reference to Figure la. Most commonly, in order to compensate 

10 for a gain variation of the output screen of the image intensifier a first corrective image is 
constructed through averaging a set of raw images, said images preferably being acquired 
using a homogeneous test phantom. Alternatively, an in-air scan can be acquired for this 
purpose. An example of a suitable method for an image correction is given by: 

to (i, j,n) =Im^(i, j,n) '-^^ -— , where 

gain^^Oj) 

15 Inigain. singiestagewith raw images is a rcsuIting gain-corrccted image of an object under investigation, 
for example a patient; 

Imraw is an original non-processed raw image of the object to be reconstructed, of a three- 
dimensional scan consisting of M images (l<n<M); 

gainrawOJ) is a gain value for a pixel (i j), (l<ij^, for a gain correction image. 

20 From Figure 2a it follows that using a single gain correction for reduction of 

the unage artifact results in residual rings Rl, visible in the slice II. In spite of the fact that 
the intensity of this ring-shaped artifact is reduced with respect to the original non-processed 
image, it is still visible in the resulting gain-corrected image. 

Another known approach of a single gain correction is to gain correct at the 

25 level of unwarped images. Figure 2b shows an image artifact post first gain correction as 

known in the art in a slice 12 corresponding to a plane perpendicular to the rotational axis le 
of the image intensifier as is set forth with reference to Figure la. This gain correction is 
carried out with a single gain image, this gain image being genemted through averaging raw 
images of the gain calibration scan which are unwarped. The unwarping operation results in 

30 eliminating the image artifact due to the pincushion and magnetic deformation of the image. 
Suitable functions to carry out the unwarping operation are known per se in the art. Gain 
correction of unwarped images is given by: 
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6 



I^gain. single stage with utiw8ip«l imascs 0J> ^) =^ TTTT » WhcrC 

S«lnunwarpe(l"»JJ 



Inigain, s{ngiestagewithimwaipedimages(ij»n) is a value of a pixel (i j) of a gain-coiTected image; 
In)un>vaiped(iJ»n) is an unwarped image, of a three-dimensional scan consisting of M images 
5 (l<n<N4); 

gainunwaipedO j) is a gain image constructed for imwarped raw images. 

From Figure 2b it follows that after a single gain correction based on 
unwarped images the ring-shaped artifacts R2 remain in the reconstructed slice 12. 

Figure 3 shows a schematic view of a block diagram of an embodiment of an 

10 X-ray examination apparatus comprising the imaging system according to the invention. The 
X-ray examination apparatus 10 comprises an imaging system 15a arranged with the image 
artifact reduction means 20 in accordance with the invention. An X-ray beam 12 emanating 
from an X-ray source (not shown) is transmitted though an object 13 imder investigation, 
most commonly a patient. The X-rays are attenuated by the patient and the attenuated 

15 radiation 12 is collected by the input screen 14 of the image intensifier. An X-ray image 

collected on the input screen of the image intensifier is converted into a visible image on the 
output screen 16 of said image intensifier. By means of an optical arrangement 17 said visible 
image is imaged onto an imaging means 18, preferably a CCD camera, whereby the visible 
image is converted into an electrical analog video signal 19. According to the invention the 

20 X-ray apparatus comprises the artifact reduction means 20, said means being arranged to 
subsequently process the image using a first function for eliminating a first structured noise 
in the image and a second function for eliminating a second structured noise in the image. 

Preferably, to account for the structured noise originating from the output 
screen of the image intensifier, a first corrective image is constructed through averaging a set 

25 of raw images acquired during a three-dimensional gain calibration scan. In order to enable 
an accurate double stage gain correction for the three-dimensional scan of, for example a 
patient, tlie geometry of the three-dimensional scans of the calibration object and the patient 
has to be substantially the same.The resulting first gain correction data for raw images is 
given by: 

1 ^ 

30 gain niw(i, j, ) = — X S^^^'^^h j\ n) , 
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where gainraw(ij,n) is the series of raw gain images of the gain calibration scan acquired 
during a three-dimensional scan; 

gainrawOJO is the first gain correction data constructed for a set of N gain calibration raw 
images comprising a z xy number of pixels. 
5 Thus, this operation involves averaging of a set of raw images. Raw image 

pixels have a one-to-one relation with corresponding elements of the CCD camera, 
independent of the image number. The elements have a fixed relation with spots on the 
output screen of the image intensifier, independent of the orientation of the image intensifier. 
However, due to pincushion deformation, a certain CCD pixel does not have a fixed relation 
10 with a spot on the input screen of the image intensifier. Therefore, high frequency structured 
noise of the input screen is averaged out when constructing gaii^wCU )- So gainnwCiJ) 
describes gain effects of the output screen of the image intensifier, the influence of the input 
screen being averaged out. 

A further method for suppressing the high-fiequency part of the structured 
15 noise of the input screen is to defocus tiie image intensifier. It may occur that a projection of 
the output screen of the CCD camera moves slightly during an acquisition of a three- 
dimensional scan, it being caused by inertial and gravitational forces which deform the 
camera assembly; such a movement can be calibrated through analyzing the drift of the 
structured noise, as observed with the CCD during the gain calibration scan. An example of a 
20 suitable drift calibration comprises a parameterization of the movement of the projection of 
the output screen onto the CCD camera in terms of a x-pixel shift, a y-pixel shift and a 
rotation of image n of the scan with respect to the image 1 of the scan. This parameterization 
is then accounted for during a construction of the first gain correction data and is applied to 
the raw images of the patient scan. 
25 In order to account for the structured noise of the input screen a second 

connective image for a set of unwarped images of the gain calibration scan is constructed. The 
corresponding gain function is given by: 
1 ^ 

gain uflwarped(i, j,) = — y\gammmarped(i^jyn) , where 

g£ununwaiped(i J^n) is the series of gain images of the gain calibration scan acquired during a 
30 three-dimensional scan, said images being corrected with the first corrective image 
gainraw(ij>) and subsequently unwarped. 

Unwarped image pixels have a fixed relation with a position on the pincushion 
grid used for a warp calibration, which is used for purposes of unwarping. Therefore, 
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unwarped image pixels have a substantially fixed relation with spots on the input screen of 
the image intensifier, independent of its orientation. Due to the magnetic image deformation, 
a certain unwarped image pixel does not have a fixed relation with a spot on the output 
screen of the image intensifier. Therefore, the structured noise of the ou^ut screen is 
5 averaged out when constructing gainunwaiped(ij)- So, gainunwaiped(ij) describes gain effects of 
the input screen of the image intensifier, not of tfie output screen. The effect of the warping is 
used to separate the structured noise of the input screen and the output screen of the image 
intensifier. 

In order to carry-out a double-stage gain correction of the images of, for 
10 example a patient scan, the first corrective image gain taw(ij»n) and the second corrective 
image gainunwarped(ij»n), as discussed above, are applied to the three-dimensional scan of the 
patient. This is discussed in detail with reference to Figure 4. The resulting image corrected 
by means of the double-stage gain correction is made available by the artifact reduction 
means 20 to a suitable visualization means, for example to a computer monitor 30 for 
15 inspection purposes. An example of a resulting reconstructed image 13 is shown in Figure 5. 
It is seen that when compared to the images shown in Figures 2a and 2b the ring-shaped 
artifact is substantially removed due to the double stage gain calibration according to the 
invention. 

The artifact reduction means 20 according to the invention preferably 
20 comprises a first stage correction means 21 arranged to eliminate the structured noise of the 
output screen using the first corrective image. Preferably, the first corrective image 21a is 
pre-calculated and is stored in a suitable memory unit of a computer (not shown). A raw 
image of the patient is first processed with the first corrective image 21a. The thus obtained 
gain-corrected image is forwarded to an image deformation correction means 23, where a 
25 suitable unwarping function 23a is being pre-stored. The resulting gain-corrected unwarped 
images are then made available to the second stage gain correction means 25, where a second 
corrective image 25a is applied to the images resulting in a final set of images with a 
substantially reduced ring-shaped artifact. The final set of images is made available to the 
image reconstruction means 26 arranged for further processing of the final set of images, the 
30 result thereof being, for example, visualized on a computer monitor 30 for inspection 
purposes. 

Figure 4 shows schematically a flowchart 40 of an operation of the computer 
program according to the invention. Prior to the double-stage calibration step C, a three- 
dimensional gain calibration scan 41 and a three-dimensional patient scan 51 are acquired. 
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The three-dimensional gain calibration scan 41 is used in a preparatory step G to obtain the 
first corrective image gainraw(ij) 43 and the second corrective image gainunwarped(ij) 45, 
respectively, according to the procedure discussed above. Preferably, the first corrective 
image 43 and the second corrective image 45 are pre-stored in respective look-up tables 52, 
5 54 and are addressed by the computer program during its execution. Thus, at step 53 the 
double-stage gain correction algorithm carries out a gain correction of the raw images of the 
patient by using the first corrective image 52, yielding the gain-corrected image of the patient 
according to the equation: 

■^Ssain„^(i,j) 

Inigain. raw J > = Im^ (1, j , u) '-^^ -— ,where 

gain^^(i,j) 

10 Imgain, raw(ij,n) is a gain-corrected image of the patient; 
hn raw(ijjn) is a raw image of the patient; 

gainraw(i J) is a first corrective fiinction constructed for a set of raw images of a gain 
calibration scan. 

Additionally, it is possible to correct for a drift of the projection of the output 
15 screen of the image intensifier of the imaging means by applying a drift correction 53b 
prestored in a look-up table 53a. 

Subsequently the gain-corrected images Imgaiii,raw(iJ,n) are subjected to an 
unwarping operation at step 55 by means of a suitable known per se in the art imwarping 
fiinction 23a stored in an image deformation correction means 23 (see Figure 3), yielding 
20 unwarped gain-corrected images of the calibration object Imgain, unwarpedCi J»n). At this stage the 
structured noise of the output screen of the image intensifier is substantially removed fi-om 
the raw images of the patient, said images being also corrected for the pincushion 
deformation caused by external magnetic fields. 

Finally, at step 57 the g^in-corrected unwarped images Imgain. unwaiped(io»i^) ^ure 
25 corrected with the second corrective image gmnunwoipedOJX preferably pre-stored in a look-up 
table 54, this being given by the equation: 

•^ZS^^ unwarped OJ) 

Imgain. doublestageOJ»Il) = ^"1 gain. unwarped 0 J» n) ^ — , Where 

g^I^ unwarped U» 

Imgain, doublestage (ij,n) is a double-stage gain corrected image of the patient; 
Imgain, imwaiped(ij,n) is an unwarped gain corrected image of the patient; 



wo 2005/006257 



PCT/IB2004/051135 



10 

gainunwaiped(i J) is the second corrective image constructed for a set of unwarped gain- 
corrected imaged of a gain calibration scan. 

The resulting image is made available by the computer program for the further 
processing at step 59, preferably it is stored in a suitable memory unit of a computer. 



